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Abstract 

A number of rice resistance genes, called Xa genes, have been identified that confer resistance against various 
strains Xanthomonas oryzae pv. oryzae (Xoo), the causal agent of bacterial blight. An understanding of pathotype 
diversity within the target pathogen population is required for identifying the Xa genes that are to be deployed for 
development of resistant rice cultivars. Among 1024 isolates of Xoo collected from 20 different states of India, 11 
major pathotypes were distinguished based on their reaction towards ten Xa genes {Xa1, Xa3, Xa4, xa5, Xa7, xa8, 
Xa10, Xa11, xa13, Xa21). Isolates belonging to pathotype III showing incompatible interaction towards xa8, xa13 and 
Xa21 and compatible interaction towards the rest of Xa genes formed the most frequent (41%) and widely distributed 
pathotype. The vast majority of the assayed Xoo isolates were incompatible with one or more Xa genes. 
Exceptionally, the isolates of pathotype XI were virulent on all Xa genes, but have restricted distribution. Considering 
the individual R-genes, Xa21 appeared as the most broadly effective, conferring resistance against 88 % of the 
isolates, followed in decreasing order by xa13 (84 %), xa8 (64 %), xa5 (30 %), Xa7 (17 %) and Xa4 (14 %). Fifty 
isolates representing all the eleven pathotypes were analyzed by southern hybridization to determine their genetic 
relatedness using the \S1112 repeat element of Xoo. Isolates belonging to pathotype XI were the most divergent. 
The results suggest that one RFLP haplotype that is widely distributed all over India and is represented in strains 
from five different pathotypes might be an ancestral haplotype. A rice line with xa5, xa13 and Xa21 resistance genes 
is resistant to all strains, including those belonging to pathotype XI. This three gene combination appears to be the 
most suitable Xa gene combination to be deployed in Indian rice cultivars. 
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Introduction 

Xanthomonas oryzae pv. oryzae (Xoo) is the causal agent of 
bacterial blight, a serious disease of rice. Bacterial blight is 
endemic and causes serious yield losses for the rice crop 
grown in irrigated, low land areas across Asia. Host plant 
resistance is the most effective way of managing yield losses 
due to the disease as chemical control is not effective [1]. 
Almost thirty different rice genes (called Xa genes) that confer 
resistance against various races and pathotypes of Xoo have 
been identified [2]. Many of these resistance genes have been 
tagged with closely linked molecular markers and are being 
used in marker assisted selection [3-6]. These resistance 
genes display specificity with regard to their effectiveness 
against different pathogen races. Therefore knowledge of the 



pathotype diversity in the target pathogen population is 
essential for making an informed choice of resistance genes 
that are to be used in a breeding program. 

In India, bacterial blight occurs in a large number of states 
with yield losses going up to 60-80% in severe infections. 
Pathotype studies conducted in India under the All India 
Coordinated Rice Improvement Project (AlCRIP) during the 
1970s and 1980s indicated that pathotypes la and lb were 
distributed all over the country [7,8]. These two pathotypes 
exhibit similar reactions on differential cultivars and are 
characterized by incompatibility with rice varieties BJ1 (xa5 and 
xa13) and DV85 (xa5, Xa7 and xa24). They are distinguished 
from each other on the basis of reaction pattern on rice cultivar 
IR20 which carries the Xa4 resistance gene; pathotype la 
strains are incompatible on IR20 while pathotype lb strains are 
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compatible. \S1112, a repeat element native to Xoo, is a good 
probe for genotyping different strains of the pathogen as it is 
present in multiple copies and reveals substantial inter-strain 
variability [9-12]. Yashitola at al. (1997) performed pathotype 
analysis and DNA fingerprinting studies with \S1112 on a set of 
Indian isolates of Xoo that were collected between 1991 to 
1995 [13]. They observed that the vast majority of strains 
(60/67) analyzed were incompatible with BJ1 and DV85 but 
were compatible with IR20. This suggested that all of these 
strains belonged to pathotype lb and none were like pathotype 
la. RFLP haplotyping demonstrated that strains belonging to 
this pathotype clustered together in a dendrogram as a lineage 
comprised of closely related strains. The BX01 Xoo strain was 
considered as the type strain for this lineage. The remaining 
strains (7/67) were found to be compatible with BJ1 and DV85 
(therefore they were neither pathotype la nor pathotype lb) and 
they had very diverse haplotypes. The BX08 strain was 
considered as an example of this group of strains. Shanti et al. 
(2001) inoculated a collection of Xoo strains from Eastern India 
on a set of near-isogenic rice lines containing any one of 
several Xa genes in the IR24 rice varietal background [14]. 
They identified several Xoo strains that are compatible with the 
xa13 disease resistance gene and suggested that the 
combination of Xa4, xa5 and Xa21, when pyramided, would be 
effective against strains from Eastern India. Lore et al. (2011) 
inoculated 224 Xoo strains collected between 1999-2006 from 
the North western Indian state of Punjab on near-isogenic rice 
lines containing individual Xa genes and found that none of 
these strains were compatible with the xa13 disease resistance 
gene [15]. Approximately 93% of the isolates were 
incompatible with Xa21 while 95% of the isolates were found to 
be compatible with the xa5 resistance gene. In an earlier study 
of Xoo strains from the Punjab, several isolates were found to 
be compatible with xa13 [16]. 

All of these above mentioned studies were done with either a 
limited number of strains or were from one particular 
geographic area of India. In order to obtain a more 
comprehensive picture of the pathotype diversity of Xoo in 
India, we have collected 1024 isolates from twenty different 
states. These strains were inoculated on a set of near-isogenic 
lines containing any one of nine Xa genes which have been 
used in earlier studies to assess pathotype diversity of Indian 
strains of Xoo (14,15). Besides these lines, the IRS rice variety 
which carries the Xa11 disease resistance gene was also 
included in the study. Thus, the effectiveness of ten different 
Xa genes was assessed and eleven different pathotypes were 
identified in this study. Representative strains from each of the 
pathotypes were subjected to DNA fingerprinting using the 
\S1112 probe and the relationship between RFLP haplotype 
and pathotype was analysed. The three gene combination of 
Xa21, xa13 and xa5 appears to be appropriate for development 
of bacterial blight resistant rice varieties in India. 

Materials and Methods 

Bacterial collection and maintenance 

Infected leaf samples were collected during 2004 - 2012, 
from different rice growing states of India representing various 



geographical locations (Figure 1). Although a majority of the 
states were sampled, because of constraints arising from the 
fact that a vast region was involved, certain states were 
sampled sparsely or were not sampled at all. All the samples 
were collected from private rice fields with owner's permission. 
Samples were placed in auto seal plastic packets with silica gel 
and stored in a refrigerator until isolation of the pathogen. 
Single colonies of cultures isolated from disease samples were 
picked up from Peptone Sucrose Agar (PSA) plates, 
maintained in liquid medium at 4°C for routine work and in 20% 
glycerol at -80°C for long term storage. 

Pathotype analysis 

The seeds of bacterial blight differential rice lines (IRBB-1, 3, 
4, 5, 7, 8, 10, 13, 21, IR8), the susceptible check line (IR24) 
and 7 different gene pyramid lines (IRBB-52, 54, 55, 57, 58, 59, 
and 60) were provided by the International Rice Research 
Institute (IRRI), Philippines. The plants were grown in plastic 
trays (55 x 40 x 15 cm) in a greenhouse. The trays were filled 
with a mixture of soil and farmyard manure at a ratio of 3:1. N- 
P-K were supplied to the plants at the rate of 100-75-74 kg/ha 
as the basal dose in the form of urea, superphosphate, and 
muriate of potash. Sowing of the differentials was done at 20 
days intervals to get 2 stages of plants (60 and 40 days old) at 
the time of inoculation. The trays were irrigated every day. 
Adequate plant protection measures were taken to ensure 
healthy and vigorous growth of the plants. Plants were clip- 
inoculated with bacterial suspensions of 10** cfu/ml [17]. Four 
leaves per plant were inoculated for each isolate-cultivar 
combination for 60 and 40 days old plants. 

Disease observations were taken 14 days after inoculation 
by measuring lesion length. Lesion lengths <5 cm were 
considered as resistant, 5-10 cm were considered as 
moderately resistant and >10 cm were considered as 
susceptible. Pathotype grouping was done based on the 
reaction pattern onto the differentials. 

Genotype analysis 

A total of 50 Xoo isolates representing all pathotypes were 
used for genomic DNA isolation and fingerprinting. Two 
previously genotyped strains from the study of Yashitola et al. 
(1997), BX01 and BX08, were also included in the present 
analysis for comparison. 

DNA isolation and Southern hybridization 

Xoo cultures were grown for 24 h in peptone sucrose 
medium on a rotary shaker at 28°C. Genomic DNA was 
isolated using the phenol-chloroform method. Genomic DNA 
was digested using restriction endonucleases EcoRI (New 
England Biolabs Inc., Ipswich, MA, USA) for 1 h at 37°C. One 
microgram of completely digested DNA from each strain was 
separated by electrophoresis on 0.7% agarose gels, 
denatured, neutralized, and vacuum-transferred to Hybond N+ 
(GE Healthcare Bio-Sciences, Uppsala, Sweden) membranes 
according to the procedure given by Sambrook et al [18]. Blots 
were pre-hybridized in a solution of 0.5 M sodium phosphate 
(pH 7.2), 7% sodium dodecyl sulfate (SDS), 1 % bovine serum 
albumin, and 1 mM EDTA for 3 h at 65°C. Probes were labelled 
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Figure 1. Statewise distribution of X. oryzae pv. oryzae strains collected from India. X. oryzae pv. oryzae strains were 
collected from different states in India. The numbers in parentheses indicate the number of strains isolated from each state. Overall, 
the states from which more isolates were collected are those that have been more extensively surveyed. 

doi; 1 0.1 371/journal.pone.0081 996.g001 



with a^^p.dATP and hybridized for 18 h at 65°C with constant 
shaking. Blots were washed three times (20 min/wash) at 65°C 
first with 2x SSC (1x SSC is 0.15 M NaCI plus 0.015 M sodium 
citrate, pH 7.0), 0.1% SDS, and 5 mM sodium phosphate (pH 
7.0) and subsequently with 0.5x SSC, 0.1% SDS, and 3 mM 
sodium phosphate buffer (pH 7.0). Autoradiography was done 
with X-ray film. A 1 kilobase DNA ladder (New England Biolabs 
Inc.) was added to all the gels as a marker. After transfer, the 
membrane was cut, and the lane carrying the marker was 
hybridized separately. 

RFLP analysis 

Twenty nine different fragments present between 1 kb and 5 
kb were used for genotyping. The haplotype (fragmenting 



pattern) of each strain obtained by Southern hybridization was 
compared with all the other strains, and presence or absence 
of a particular fragment was recorded as 1 and 0, respectively. 
Similarity matrix was prepared using the dice coefficient option 
and dendrograms were prepared using the UPGMA 
(unweighted pair group method of averages) option of the 
software FREETREE [19]. The confidence limits of 
dendrograms were determined by bootstrap analysis with 
2,000 replications using the same program. The bootstrap 
values are expressed as a percentage of these 2,000 
replications. The dendrogram was drawn using the software 
TreeView [20]. 

Using the same data, a minimum spanning tree (MST) was 
also constructed for phylogenetic analysis using the MST gold 
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Table 1. Pathotyping of Indian isolates of X. oryzae pv. 
oryzae. 



Host 

Differential^ 



X. oryzae pv. oryzae pathotypes 







II 


III 


IV 


V 




VI 


VI 




VIII 


IX 


X 




XI 


IRBB1 (Xa1) 


S 


s 


s 


s 


s 




s 


s 




s 


s 


R, 


M 


s 


IRBB3 (Xa3) 


S 


s 


s 


s 


s 




s 


s 




s 


s 


R 




s 


IRBB4 (Xa4) 


s 


s 


s 


s 


R, 


M 


M 


R, 


M 


R, M 


s 


R 




s 


IRBB5 (xa5) 


s 


s 


s 


s 


R, 


M 


R 


R 




R 


R 


R, 


M 


s 


IRBB7 (Xa7) 


s 


s 


s 


s 


S 




S 


S 




R 


R 


R, 


M 


s 


IRBB8 (xa8) 


s 


s 


R, M 


R, M 


R, 


M 


R 


R, 


M 


S 


S 


R 




s 


IRBB10 (XalO) 


s 


s 


S 


S 


S 




S 


S 




S 


S 


R, 


M 


s 


IRBB13 {xa13) 


R 


R 


R 


R 


R 




S 


S 




S 


S 


R, 


M 


s 


IRBB21 (Xa2*) 


R, M 


S 


R, M 


S 


R 




M 


R, 


M 


R 


R, M 


R 




s 


IRS (Xall) 


S 


8 


S 


S 


S 




R 


S 




S 


S 


R, 


M 


s 


IR24 


S 


S 


S 


S 


S 




S 


S 




S 


S 


R, 


M 


s 



^ Based on mean lesion lengths, the nature of responses were classified into 
susceptible (S) {above 10 cm), moderately resistant (MR) (between 5-10 cm) and 
resistant (R) (up to 5 cm). The data presented are from experiments with 60 day 
old plants; similar results were obtained from inoculation of 40 day old plants. R, M 
is indicated when some isolates in the pathotype give a resistant reaction while 
other isolates give a moderately resistant reaction 

^ The Xa resistance gene in each host differential is given in parenthesis. IRS is a 
cultivar with the Xa11 resistance gene and 1R24 is the susceptible control, 
doi: 10.1371/journal.pone.0081996.t001 

software programme available at 

www.bellinghamresearchinstitute.com /software [21]. The pair 
wise distances between strains were calculated using the 
equidistant method option in the software programme [21]. A 
1000 bootstrap iterations were done on 500 unique MSTs and 
the MST with highest average bootstrap percentage was taken 
as the representative MST. The MST was visualized using the 
software GVEdit for Graphviz version 1.01 [22]. Using the 
same software, a minimum spanning consensus network was 
also made, using the pathotyping data in Table 1, for the 
eleven pathotypes identified in this study. A binary score was 
generated by taking either a resistant or a moderately resistant 
interaction as 1 and a susceptible interaction as 0. 

Results 

Pathotyping of X. oryzae pv. oryzae strains 

A total of 1024 isolates of Xoo were collected from 20 
different states in India (Figure 1) during 2004 to 2012. All of 
these isolates were pathotyped using 9 NILs harbouring Xoo 
resistance genes Xa1, Xa3, Xa4, xa5, Xa7, xa8, Xa10, xa13, 
Xa21 and one variety (IRS) that is reported to carry the Xall 
resistance gene. This analysis revealed 1 1 distinct pathotypes 
(Table 1; Figure 2). The pathotype III was the most frequent 
pathotype and accounted for 40.7 % of the isolates. Strains 
belonging to this pathotype were widely distributed throughout 
India with presence in 19/20 sampled states in the country 
(Table 2). Pathotype I was the second most prevalent 
pathotype comprising 20.5 % of the isolates that were 



45.0(1 

40.00 

35.00 

£. 30.00 

>* 

u 

I 25.00 

0) 20.00 
li. 

15.00 
10.00 
5.00 
0.00 



1 



I I I.I 



I II III I\ \ M VII VIII L\ X XI 

Pathotypes 

Figure 2. Frequency of X. oryzae pv. oryzae pathotypes in 
India. Eleven different pathotypes, designated as I to XI, were 
identified in this study. The x-axis indicates pathotype and the 
y-axis indicates % frequency amongst the 1024 isolates that 
were collected and analyzed. 

del: 1 0.1 371/journal.pone.OOSI 996.g002 

analyzed. This pathotype was present in 14 of the states that 
were covered in this study. As compared to these two 
pathotypes, the remaining nine pathotypes were isolated less 
frequently with each of them being present in this study sample 
at a frequency that was <10%. Pathotype X exhibited an 
incompatible interaction with all differentials as well as 
susceptible control variety (IR24) indicating that this pathotype 
consists of less virulent strains. This pathotype was present in 
11/20 sampled states and constituted 4 % of the total Xoo 
isolates. It is likely that pathotype X is a grouping of strains 
from other pathotypes that have, for unknown reasons, lost 
their virulence. Pathotype XI consisted of six isolates collected 
from Tripura in North-East India and showed a compatible 
interaction with all the differentials. Except for pathotype XI, all 
other pathotypes that are compatible with xa13 are 
incompatible with xa5 and those that are compatible with xa5 
are incompatible with xa13. Pathotypes V and X are 
incompatible with both xa5 and xa13. Strains from all 
pathotypes, except for pathotype XI, are found in the Eastern 
Indian state of Odisha (Table 2). The pathotype III is the most 
frequent pathotype in 16/20 states covered in this study. 
Interestingly, pathotypes that are compatible with xa13 were 
not isolated in half of the states that have been surveyed. 
These include the states of Andhra Pradesh, Karnataka, 
Maharashtra, Gujarat, Himachal Pradesh, Uttarakhand, 
Jharkand, Kerala, Manipur and Madhya Pradesh. However, the 
sample size was limited in some of these states. Although a 
large number of isolates were collected from the states of 
Haryana and Punjab, only a few isolates (1 out of 72 isolates 
and 3 out of 130 isolates respectively) were compatible with 
xa13. In comparison, many of the isolates that are compatible 
with xa13 were found in the states of Eastern India, particularly 
in Odisha and Chhattisgarh (54 out of 152 and 25 out of 193, 
respectively). 

Considering the individual R-genes, Xa27 appeared as the 
most broadly effective, conferring resistance against 88 % of 
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Table 2. Distribution of X. oryzae pv. oryzae patliotypes in different states of India. 



state' 



Pathotypes 



Total isolates 





1 


II 


III 


IV 


V 


VI 


VII 


VIII 


IX 


X 


XI 




Andhra Pradesh 


14 (16.5) 


4(4.7) 


56 (65.9) 


3 (3.5) 


5 (5.9) 










3 (3.5) 




85 


Assam 






1 (14.3) 




4 (57.1) 






1 (14.3) 


1 (14.3) 






7 


Bihar 


10 (18.2) 


9 (16.4) 


12 (21.8) 


1 (1.8) 


2 (3.6) 








21 (38.2) 






55 


Chhattisgarh 


32 (16.6) 


1 (0.5) 


96 (49.7) 


4(2.1) 


23 (11.9) 




17(8.8) 




8(4.1) 


12(6.2) 




193 


Gujarat 






12 (70.6) 




5 (29.4) 














17 


Himachal Pradesh 






4 (57.1) 


1 (14.3) 


2 (28.6) 














7 


Haryana 


12 (16.7) 


3 (4.2) 


32 (44.4) 


13 (18.1) 


8 (11.1) 






1 (1.4) 




3 (4.2) 




72 


Jharkhand 






6 (100) 


















6 


Karnataka 


1 (6.7) 




10 (66.7) 














4 (26.7) 




15 


Kerala 


1 (12.5) 


2 (25.0) 


4 (50.0) 


1 (12.5) 
















8 


Manipur 






2 (66.7) 




1 (33.3) 














3 


Maharashtra 


17 (38.6) 


5 (11.4) 


18 (40.9) 


3 (6.8) 


1 (2.3) 














44 


Madhya Pradesh 




1 (14.3) 


5 (71.4) 


1 (14.3) 
















7 


Odisha 


26 (17.1) 


20 (13.2) 


34 (22.4) 


10 (6.6) 


4 (2.6) 


5 (3.3) 


22 (14.5) 


10 (6.6) 


17 (11.2) 


4 (2.6) 




152 


Punjab 


36 (27.7) 


9 (6.9) 


59 (45.4) 


3 (2.3) 


16 (12.3) 




1 (0.8) 




2(1.5) 


4(3.1) 




130 


Tamil Nadu 


2 (25.0) 




4 (50.0) 












2 (25.0) 






8 


Tripura 


1 (7.1) 














3 (21.4) 


4 (28.6) 




6 (42.9) 


14 


Uttar Pradesh 


24 (25.0) 


3(3.1) 


36 (37.5) 


4 (4.2) 


4 (4.2) 




14(14.6) 


1 (1.0) 


10 (10.4) 






96 


Uttarakhand 


27 (39.7) 


5 (7.4) 


16 (23.5) 


8 (11.8) 


1 (1.5) 










11 (16.2) 




68 


West Bengal 


7(18.9) 


1 (2.7) 


10 (27.0) 


3(8.1) 








2 (5.4) 


14 (37.8) 






37 


1024 



^ The number of isolates of each pathotype is provided along with the total number of isolates collected from the state . The numbers in parenthesis indicate % frequency, 
doi: 10.1371/journal.pone.0081996.t002 



the isolates, followed In decreasing order by xa13 (84 %), xaS 
(64 %), xa5 (30 %), Xa7 (17 %) and Xa4 (14 %) (Figure 3). The 
R-genes Xa1, Xa3 and Xa10 were not effective against any of 
the pathotypes, except for those strains belonging to pathotype 
X. 

The reaction of the isolates against a set of gene pyramid 
lines carrying several combinations of 'Xa' genes was 
assessed (Table 3). The IRBB52 line which harbours both Xa4 
and Xa21 genes exhibited a moderate level of resistance 
against pathotype II and a resistant reaction against pathotype 
VI. The NILS that carry either one of these individual resistance 
genes exhibited a susceptible reaction with pathotype II and a 
moderate level of resistance against pathotype VI. In particular, 
the reaction with pathotype II indicates that the presence of 
both of the 'defeated' genes leads to an enhanced level of 
resistance, a phenomenon that has been previously described 
and termed as quantitative complementation (QC) [3,4,23]. 
Similarly QC was also observed in IRBB 52, IRBB54 (xa5 + 
Xa21), IRBB55 {xa13 + Xa21) and IRBB 58 (Xa4 + xa13 + 
Xa21) lines which exhibited a moderate level of resistance 
against pathotype XI as compared to the susceptible reaction 
in lines carrying the respective single resistance genes. In 
another example of QC, IRBB54 exhibited a moderate level of 
resistance towards pathotype IV as compared to the 
susceptible interaction exhibited by lines carrying either xa5 or 
Xa21. In an important observation, the IRBB59 line (xa5 + xa13 
+ Xa21) is found to exhibit resistance to pathotype XI even 
though the lines carrying each of the individual genes are 
susceptible. Interestingly, the IRBB55 line exhibits a 



^70 
^60 



Resistant 
I I Moderately resistant 




Xal Xa3 Xa4 xaS Xa7 sa8 XalO 



13 Xa2l IKS IK24 



Resistance genes 

Figure 3. Effectiveness of Xa genes against Indian X. 
oryzae pv. oryzae isolates. The x-axis indicates resistance 
gene and the y-axis indicates the frequency of isolates against 
which this gene confers resistance. Inoculation of rice leaves 
with isolates was done as described in Methods and lesion 
lengths were measured 14 days after inoculation. Lesion 
lengths <5 cm are considered as resistant, 5 to 10 cm are 
considered as moderately resistant and >10 cm are considered 
as susceptible. 

doi: 10.1371/journal. pone. 0081 996.g003 
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Table 3. Reactions of X. oryzae pv. oryzae pathotypes on 
Xa gene pyramid lines. 





Gene Pyramids^ 


Xoo pathotypes ^ 
















1 II III IV V 


VI 


VII 


VIII IX 


X 


XI 


IRBB52(Xa4 + Xa21) 


R M° R S R 


r'= 


R 


R 


R 


R 


M= 


JRBB54 {xaS + Xa21) 


R S R M° R 


R 


R 


R 


R 


R 


M= 


IRBB55 (xa13 + Xa21) 


R R R R R 


s" 


M 


R 


M 


R 


M= 


IRBB58 (Xa4 + xa13 + Xa21) 


R R R R R 


M 


R 


R 


M 


R 


M° 


IRBB59 {xaS + xa13 + Xa21) 


R R R R R 


R 


R 


R 


R 


R 


R= 


IRBB60 (Xa4 + xaS + xai3 + Xa21) 


R R R R R 


R 


R 


R 


R 


R 


R° 



^ Based on mean lesion length, the nature of responses were classified into 
susceptible (S) (above 10 cm), moderately resistant (M) (between 5-10 cm) and 
resistant (R) {up to 5 cm). 



° The Xa genes that are present in the gene pyramid line is indicated in 
parentheses. 

^ The gene pyramid lines exhibit more resistance than the lines having the single 
'R' genes. 

^ The gene pyramid line exhibits less resistance than the lines having the single 'R' 
genes. 

doi: 10.1371/journal.pone.0081996.t0C3 

susceptible response towards pathotype VI while the line 
carrying Xa21 alone exhibits a moderate level of resistance 
against this pathotype. This suggests that, in the presence of 
xa13, there is a reduction in the level of resistance that is 
provided by Xa21 against pathotype VI. 

In order to determine the shortest path by which the eleven 
pathotypes could be derived from each other, we have 
performed a minimum spanning network analysis (Figure SI) 
using the pathotyping data from Table 1. We have added, to 
the different edges (connecting lines) in this network, a 
hypothetical directionality by assuming that in most cases the 
evolution of a pathotype would have involved acquisition of the 
ability to break down a host resistance gene. Pathotype V is 
incompatible with the xa5 and the xa13 resistance genes. 
Therefore, pathotype V could be an ancestral pathotype from 
which strains that are compatible with xa5 or xa13 could have 
evolved by acquisition of the ability to breakdown these 
resistance genes. The hypothesized change from pathotype V 
to pathotype III might have involved acquisition of the ability to 
break down the xa5 resistance gene as well as Xa4. Pathotype 
II might have arisen from pathotype I through acquisition of the 
ability to break down the Xa21 resistance gene or from 
pathotype IV through acquisition of the ability to break down 
the xaS resistance gene. Pathotype XI may have arisen from 
pathotype II through acquisition of the ability to breakdown the 
xa13 resistance gene. The hypothesized change from 
pathotype V to pathotype VII might have involved acquisition of 
the ability to break down the xa13 resistance gene. The 
hypothesized change from pathotype VII to VI would have 
resulted in susceptibility to the Xa11 resistance gene and the 
change from pathotype VII to VIII would have required 
acquisition of the ability to break down the xa8 resistance gene. 
In comparison to pathotype VII, pathotype VIII is also avirulent 
on the Xa7 resistance gene. Pathotype IX may have arisen 



from pathotype VIII through acquisition of the ability to 
breakdown the Xa4 resistance gene. 

RFLP analysis Xanthomonas oryzae pv. oryzae 
strains 

RFLP based genotyping of 50 isolates belonging to 11 
different pathotypes was done by scoring 29 bands obtained by 
using \S1112, an Xoo insertion sequence element, as the DNA 
fingerprinting probe. Five strains each of pathotypes V, VIII, IX 
and X; four strains each of pathotypes I, III, IV and VI; six 
strains each of pathotypes II and VII and two strains of 
pathotype XI were analyzed (Table 4). Two other strains, 
namely BX01 and BX08, which were pathotyped and 
genotyped in an earlier study [13] were also included as 
reference strains. The RFLP and phylogenetic analysis using 
UPGMA revealed that 20 haplotypes were present amongst the 
50+2 isolates that were studied (Figures 4 and Figure S2). The 
most frequent RFLP haplotype consisted of 13 isolates 
including the following: 4 strains belonging to pathotype III, 
each of which was isolated from four different states (Andhra 
Pradesh, Chhattisgarh, Punjab and West Bengal), 4 strains 
belonging to pathotype V isolated from different states (Andhra 
Pradesh, Bihar, Maharashtra and Punjab) as well as two 
strains each of pathotypes I and IV and one strain from 
pathotype X. This indicates that pathotypes III as well as V 
consisted of genetically closely related strains which had 
dispersed to several widely distributed locations in India. The 
BX01 strain, which was included as a reference strain, is found 
to have the same RFLP haplotype as these 13 strains. The 
study of Yashitola et al. (1997) had indicated that strains with 
the same RFLP haplotype as BX01 were widely distributed in 
India. The BX01 strain belongs to pathotype III. Isolates of 
pathotype VI were isolated from a specific region in India (the 
state of Odisha) and all the four genotyped isolates of this 
pathotype have a single haplotype. Isolates of pathotype IX 
that were isolated from Odisha (East), Uttar Pradesh (North) 
and Tamil Nadu (South) clustered together in the dendrogram. 
In contrast isolates belonging to pathotypes I, II and X were 
found to be composed of genetically diverse strains as they 
grouped into different clusters. Genotyping of two strains 
belonging to pathotype XI indicated that they belong to a 
haplotype which is quite divergent from the rest of the strains. 
Also, the BX08 strain which had been previously found to be a 
very diverse strain [13], is also found to be an outlier in this 
study. 

Additional phylogenetic analysis was performed, using the 
same RFLP fragment data, by constructing a minimum 
spanning tree which can help in hypothesizing on ancestral or 
emergent RFLP haplotypes (Figure 5). The RFLP haplotype 
representing 13 strains (14 including BX01) occupies an 
approximately central position in the tree. As indicated above, 
this is a widely distributed RFLP haplotype in India and 
includes strains from five different pathotypes. These 
observations, taken together, suggest that this might be an 
ancestral RFLP haplotype. Interestingly, this RFLP haplotype 
includes 4 strains of pathotype V which is hypothesized to be 
an ancestral pathotype. An RFLP haplotype that is shared by 7 
strains (four of pathotype VI and one each of pathotypes VII, 
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Table 4. List of representative X. oryzae pv. oryzae strains 
that were used for Restriction Fragment Length 
Polymorphism Analysis. 



Year of 



SI. No. 


Isolate 


Location 


Rice cultivar^ 


isolation 


Pathotype 


1 


1X0221 


Andhra Pradesh 


MTU1010 


2005 


1 


2 


1X0278 


Uttar Pradesh 


Arize 6444 


2005 


1 


3 


1X0812 


Andhra Pradesh 


Arize 6516 


2007 


1 


4 


1X0675 


Punjab 


Pusa 44 


2007 


1 


5 


1X089 


Uttarakhand 


NA 


2004 




6 


1X0884 


Maharashtra 


MTU 7029 


2008 




7 


1X0685 


Kerala 


NA 


2006 


II 


8 


1X0725 


Maharashtra 


H93025 


2007 


II 


9 


1X041 1 


Chhattisgarh 


Mahamaya 


2006 


II 


10 


1X0651 


Punjab 


Pusa Basmati 1 


2006 


II 


11 


1X0639 


Punjab 


Pusa 44 


2006 


III 


12 


1X0281 


Andhra Pradesh 


MTU1001 


2005 


III 


13 


1X074 


Chhattisgarh 


Arize 6444 


2004 


III 


14 


1X0189 


West Bengal 


H93024 


2005 


III 


15 


1X092 


Punjab 


PR118 


2005 


IV 


16 


1X0367 


Uttarakhand 


NA 


2006 


IV 


17 


1X035 


Chhattisgarh 


IRBB 5 


2004 


IV 


18 


1X098 


Himanchal 
Pradesh 


NA 


2005 


IV 


19 


1X099 


Andhra Pradesh 


BPT 5204 


2005 


V 


20 


1X0518 


Gujarat 


Arize 6444 


2006 


V 


21 


1X0414 


Bihar 


Arize 6129 


2006 


V 


22 


1X0134 


Maharashtra 


NA 


2005 


V 


23 


IXO220 


Punjab 


Arize 61 1 1 


2005 


V 


24 


IXO620 


Odisha 


MTU 7029 


2005 


VI 


25 


1X0644 


Odisha 


NA 


2005 


VI 


26 


1X0621 


Odisha 


Surendra 


2005 


VI 


27 


IXO630 


Odisha 


NA 


2005 


VI 


28 


1X0222 


Odisha 


NA 


2005 


VII 


29 


1X0599 


Uttar Pradesh 


Arize 6444 


2005 


VII 


30 


1X0597 


Uttar Pradesh 


Arize 6444 


2005 


VII 


31 


1X0631 


Odisha 


NA 


2005 


VII 


32 


IXO390 


Odisha 


NA 


2005 


VII 


33 


1X0792 


Chhattisgarh 


MTU 7029 


2007 


VII 


34 


IXO90 


Uttar Pradesh 


BPT 5204 


2004 


VIII 


35 


1X0493 


Odisha 


NA 


2006 


VIII 


36 


1X0169 


Uttar Pradesh 


Arize 6444 


2005 


VIII 


37 


1X0141 


Odisha 


Surendra 


2005 


VIII 


38 


1X0365 


Odisha 


MTU 7029 


2005 


VIII 


39 


1X093 


Uttar Pradesh 


Indrasan 


2005 


IX 


40 


IXO608 


Odisha 


MTU 7029 


2007 


IX 


41 


IXO603 


Odisha 


MTU 7029 


2007 


IX 


42 


1X0159 


Uttar Pradesh 


Roopali 


2005 


IX 


43 


1X097 


Tamil Nadu 


NA 


2005 


IX 


44 


1X053 


Karnataka 


NA 


2004 


X 


45 


1X0151 


Haryana 


HKR 47 


2005 


X 


46 


1X0842 


Chhattisgarh 


Arize 6444 


2007 


X 


47 


1X0645 


Odisha 


Arize Tej 


2005 


X 


48 


IXO704 


Punjab 


Rasi 7001 


2007 


X 


49 


1X01 088 


Tnpura 


MTU 7029 


2009 


XI 


50 


1X01 104 


Tripura 


MTU 7029 


2009 


XI 



Table 4 (continued). 



^ Rice cultivar from which the strain was isolated: NA indicates that the information 
is not available 

doi: 10.1371/joumal. pone. 0081 996.t004 

VIII and X) is connected to this hapiotype and couid be derived 
from it. Aiso, the postulated ancestral RFLP hapiotype is 
connected by an edge having a strong bootstrap value to an 
RFLP hapiotype shared by two pathotype VII strains 
suggesting that the latter might be an emergent hapiotype. A 
different RFLP hapiotype (strain 1X093.9) is connected to five 
different haplotypes. Support for this grouping is provided by 
the observation that two nodes representing three strains 
belonging to pathotype IX are connected to 1X093.9. 
Interestingly, two pathotype II strains that share an RFLP 
hapiotype are connected to 1X093.9 through an edge that has 
a high bootstrap value. As was observed in the tree made by 
UPGIVIA, the pathotype II strains appear at multiple places in 
the minimum spanning tree indicating that they are a diverse 
set of strains. Further studies on pathotype II strains using 
additional rice differentials might reveal significant pathological 
differences between them. 

Discussion 

Based on the study sample of 1024 strains, and the 
differential rice lines used in this study, eleven different 
pathotypes of this bacterium have been identified. Although a 
number of pathotypes were identified, pathotype XI isolated 
from the Eastern Indian state of Tripura was very unique in that 
it was compatible with all the tested Xa genes. Fortunately, the 
gene pyramid line carrying Xa21, xa13 and xa5 is resistant to 
this strain. In view of the possible dispersal of this pathotype to 
other locations in India, efforts should be made to deploy these 
three resistance genes in the genetic background of important 
Indian rice cultivars. 

Genotype analysis, using the \S1112 probe, carried out on a 
subset of 50 strains that are drawn from the 11 pathotypes 
(plus the reference strains BX01 and BX08) revealed 20 RFLP 
haplotypes indicating that the genetic base of the Indian Xoo 
population is fairly diverse. At the same time, the bootstrap 
values for many of the clusters were not significant. Multiple 
RFLP probes, PGR based finger printing methods and 
genotyping using single nucleotide polymorphisms (SNPs) 
might be needed to delineate the phylogenetic relationships 
between these strains in a more accurate manner. Using the 
same RFLP data, a phylogenetic tree was also constructed 
using minimum spanning tree analysis. One particular RFLP 
hapiotype shared by thirteen different isolates plus BX01 could 
be an ancestral hapiotype as it is centrally located in the tree 
and five different pathotypes have this hapiotype. 

Although it might be considered as being too simplistic, we 
have tried to use the data in Table 1 to generate a minimum 
path by which these eleven pathotypes could have arisen. As 
pathotype V strains are incompatible with both of the xa5 and 
xa13 disease resistance genes, pathotypes that are compatible 
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Figure 4. Restriction Fragment Length Polymorphism 
analysis of X. oryzae pv. oryzae strains. Genomic DNA was 
isolated from 52 X. oryzae pv. oryzae strains and southern 
analysis was performed using the \S1112 probe as described 
in Methods. These 52 strains included 50 strains from eleven 
pathotypes identified in this study as well as X. oryzae pv. 
oryzae strains BX01 and BX08 from the study of Yashitola et 
al (13). Twenty different haplotypes were identified and at least 
one representative of each haplotype is shown. The 1X0 strain 
number is indicated on each lane along with pathotype. A one 
l<b DNA ladder was added as a size marker. X. oryzae pv. 
oryzae strain BX01 is loaded in the second lane of each blot 
for comparison. Each strain was analyzed at least three times, 
and 29 bands that were consistently visible in all the three 
replicates were used for scoring. 

doi: 10.1371/joumal.pone.0081996.g004 



with either one or the other of these two resistance genes could 
have arisen by the acquisition of specific TAL (transcription 
activator-like) effectors by an ancestral pathotype V strain. The 
Xa5 gene encodes transcription factor TFIlAy and is 
presumably needed by Xoo TAL effectors to upregulate 
specific host susceptibility genes [24-26]. The xa5 disease 
resistance gene probably does not support this role of TAL 
effectors in up-regulating expression of host susceptibility 
factors. Xoo strains that can breakdown xa5 mediated 
resistance are postulated to produce TAL effector/s that can 
utilize an alternate TFIlAy gene in rice [25]. Therefore a 
pathotype (such as pathotype III) that can breakdown xa5 
mediated resistance could have arisen from an ancestral 
pathotype V like strain by acquisition of TAL effector/s that can 
utilize the alternate TFIlAy. Pathotypes I, II, III and IV differ 
from pathotype V in their ability to overcome xa5 as well as 
Xa4 mediated resistance. We speculate that there might have 
been some bridging strains that were able to overcome xa5 
mediated resistance but were unable to overcome Xa4 
mediated resistance. The previously described pathotype la 
which is reported to be incompatible with Xa4 could be such a 
pathotype. 

The xa13 disease resistance gene represents a mutation in 
the promoter of Os8N3, a host susceptibility gene that encodes 
a sugar transporter (SWEET protein) [27-29]. Xoo strains that 
overcome xa13 mediated resistance are capable of producing 
a TAL effector that can upregulate the expression of an 
alternate SWEET gene of the host [30]. It is possible that 
pathotype VII could have arisen from pathotype V by 
acquisition of such a TAL effector. Although it is speculative, 
we postulate that the first step in the evolution of Pathotypes 
VI, VII, VIII and IX (that can all breakdown xa13 mediated 
resistance) would have been the acquisition of such a TAL 
effector by an ancestral pathotype V like strain. 

Does the minimum spanning tree in Figure 5 provide any 
support for the notion that many of the current day Indian 
pathotypes of Xoo could have arisen from an ancestral 
pathotype V like strain? The RFLP haplotype that has been 
suggested above to be an ancestral haplotype includes four out 
of five pathotype V strains analyzed in this study. Moreover, 
several pathotype III, pathotype IV and pathotype I strains 
share the same RFLP haplotype. Therefore, it is possible that 
pathotypes III, IV and I could have evolved from an ancestral 
pathotype V strain. Figure 5 also suggests that RFLP 
haplotypes that represent pathotype VII could have emerged 
from this supposedly ancestral RFLP haplotype. Therefore, it is 
possible that pathotype VII strains could have arisen from an 
ancestral pathotype V like strain. Figure 5 also suggests that at 
least some pathotype VIII and IX strains could have arisen from 
pathotype VII strains. 

The suggestion that pathotype V might be an ancestral 
pathotype from which other pathotypes arose through 
acquisition of compatibility against either xa5 or xa13 provides 
an explanation for the observation that the pathotypes that are 
compatible with xa5 are incompatible with xa13 and vice versa. 
The exceptions include pathotype X strains which are 
incompatible with all the rice lines tested in this study. We 
suggest that pathotype X strains might have lost virulence as 
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Figure 5. Phylogenetic relationships of X. oryzae pv. oryzae strains analyzed using a Minimum spanning tree. A minimum 
spanning tree was drawn using the RFLP data. The digits in the edges represent percentage boot strap values obtained after 1000 
iterations. A solid bold edge indicates a bootstrap value over 90%; solid line indicates a value over 70% and a dotted line indicates a 
value below 70%. 

doi: 10.1371/joumal.pone.0081996.g005 



growth ceased in dried lesions of rice leaves. Spontaneous 
loss of virulence has been previously reported from aging 
cultures of Xoo [31]. Pathotype XI strains are compatible with 
both xa5 and xa13 disease resistance genes; one possible way 
in which this could have occurred is through acquisition of the 
ability to overcome xa13 mediated resistance by a pathotype 
(such as pathotype II) that is already compatible with xa5. An 
alternative possibility that pathotype XI strains arose through 
acquisition of the ability to overcome xa5 mediated resistance 
by a strain that is already compatible with xa13 cannot be ruled 
out and indeed does have some support from the minimum 
spanning tree analysis (Figure 5). 



A number of strains belonging to pathotypes III and V had 
the same haplotype even though they were isolated from 
widely separated locations. This suggests the possibility that 
they have been dispersed through seed. Pathotype IX strains 
cluster together in the dendrogram indicating that they are 
closely related strains. These strains were isolated from three 
geographically separated states in India; an observation that is 
again suggestive of the possibility of dispersal through seed. 
The BX01 strain of Xoo has been previously described by 
Yashitola et al. (1996) as belonging to pathotype lb which 
belongs to a widely distributed lineage of this pathogen in India. 
In that study, the same RFLP haplotype as that of the BX01 
strain was found in 15/67 strains that had been analyzed with 
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the \S1112 probe and these strains were widely distributed in 
India. The results presented here indicate that the BX01 strain 
has the same RFLP haplotype as the most prevalent RFLP 
haplotype found In this study. The current study confirms the 
results of Yashitola et al. (1996) and indicates that strains with 
the same RFLP haplotype as BX01 continue to be widely 
distributed in India. The same RFLP haplotype had also been 
found to be present in a majority of Xoo isolates collected from 
a wild rice species (Oryza nivara) growing naturally in Southern 
India and It had been suggested that this haplotype might have 
transferred from wild rice to cultivated rice [32]. 

Among the Xoo resistance genes studied, Xa21 was found to 
be most effective towards Indian Xoo strains followed by xa13 
and xa8. All three genes appear to be good candidates to be 
deployed in Indian rice cultivars. The current study also throws 
light on the suitability of different Xa gene combinations for 
deployment in India. Since the most effective 'R' genes are 
Xa21 and xa13, a combination of these two genes will be a 
natural choice. However, as per this study, a number of Indian 
Xoo strains that are compatible with xa13 are incompatible with 
xa5. Therefore, a gene combination of Xa21, xa13 and xa5 
might be an added advantage. This gene combination is also 
effective against pathotype XI. Two different research groups in 
India have pyramided the Xa21 and xa13 resistance genes into 
the popular rice varieties Pusa Basmati-1 and Samba Mahsuri 
[5,6]. In the latter work, the xa5 resistance gene was also 
incorporated along with Xa21 and xa13 to create a three gene 
pyramid line. 

The present study also shows examples of a phenomenon 
called quantitative complementation (QC) wherein the 
presence of two different 'Xa' genes provides an increased 
level of resistance as compared to either of the single 
resistance genes. In this study, as described in Table 3, 
several such examples of QC could be observed. Interestingly, 
we have also observed an example of antagonistic interaction 
between Xa genes, wherein the pyramid o1 Xa21 and xa13 was 
less resistant than the line having Xa21 alone. This 
antagonistic interaction was observed when this gene pyramid 
was inoculated with isolates belonging to pathotype VI. A 
similar example of pathotype specific antagonistic interaction 
between Xa21 and xa13 has been previously described [33]. A 
negative interaction between Xa21 and xa5 only in the genetic 
background of a specific rice cultivar has also been previously 
described [34]. These examples highlight the need to take into 
consideration the possible effects of genetic background of the 
cultivar and pathotype prevalence in the region while making 
decisions on deployment of Xa gene pyramid lines in bacterial 
blight affected rice growing regions. 

References 

1. Singh R, Das B, Ahmed K, Pal V (1980) Chemical control of bacterial 
leaf blight of rice. International Journal of Pest Management 26: 21-25. 
doi: 1 0. 1 080/09670878009414278. 

2. Liu-Nino DO, Ronald PC, Bogdanove AJ (2006) Xanthomonas oryzae 
pathovars: model pathogens of a model crop. Mol Plant Pathol 7: 
303-324. doi:10.1111/j.1364-3703.2006.00344.x. PubMed: 20507449. 

3. Yoshimura S, Yoshimura A, Iwata N, McCouch SR, Abenes ML et al. 
(1995) Tagging and combining bacterial blight resistance genes in rice 
using RAPD and RFLP markers. Molecular Breeding 1: 375-387. doi: 
10.1007/BF01248415. 



In summary, 1024 Xoo strains were collected from 20 
different states in India and subjected to pathotyping and RFLP 
analysis. This has provided interesting insights into the genetic 
and pathotypic diversity of Indian strains of Xoo. However, 
several states were either not sampled or were sampled rather 
sparsely. Also, very little information is available about the TAL 
effectors and the host susceptibility genes that are used by 
Indian strains of Xoo to break down rice resistance genes such 
as xa5, xa13, etc. Future studies should be aimed at 
addressing these issues. 

Supporting Information 

Figure S1. A Minimum spanning network of eleven 
different X. oryzae pv. oryzae pathotypes. The network was 
developed using data from Table 1 as indicated in methods. 

The direction of the hypothetical change is indicated by the 
arrow. The resistance gene against which compatibility is 
gained (+) or lost (-) during the change from one pathotype to 
another is indicated. Dotted lines are given when alternate 
edges are possible. The digits given on the dotted edges 
indicate the percentage of minimum spanning trees having that 
particular edge. 
(TIF) 

Figure S2. Dendrogram of 52 Indian strains of X. oryzae 
pv. oryzae derived from restriction fragment length 
polymorphism analysis using the \S1112 repeat element. 

The dendrogram was constructed and bootstrap values 
calculated as described in methods. The digits in the nodes 
represent percent boot strap values after 2000 iterations. The 
1X0 number along with the pathotype is indicated for each 
strain. The scale bar represents genetic divergence. The BX01 
and BX08 strains previously described by Yashitola et al. 1997 
(13) were included for comparison. 
(TIF) 

Author Contributions 

Conceived and designed the experiments: DIVI MGA RVS. 
Performed the experiments: DM MRV MGA KK YR. Analyzed 
the data: DM MRV MGA RVS. Contributed reagents/materials/ 
analysis tools: DM MGA RVS. Wrote the manuscript: MGA DM 
RVS. 



4. Huang N, Angeles E, Domingo J, Magpantay G, Singh S et al. (1997) 
Pyramiding of bacterial blight resistance genes in rice: marker-assisted 
selection using RFLP and PCR. Theoretical and Applied Genetics 95: 
31 3-320. doi:1 0.1 007/S001220050565. 

5. Joseph M, Gopalakrishnan S, Sharma R, Singh V, Singh A et al. (2004) 
Combining bacterial blight resistance and Basmati quality 
characteristics by phenotypic and molecular marker-assisted selection 
in rice. Molecular Breeding 13: 377-387. doi:10.1023/B:MOLB. 
0000034093.63593.4c. 



PLOS ONE I www.plosone.org 



10 



November 2013 | Volume 8 | Issue 1 1 | e81996 



Diversity of indian Xanthomonas oryzae Isolates 



6. Sundaram RM, Vishnupriya MR, Biradar SK, Laha GS, Reddy GA et al. 
(2008) Marker assisted introgression of bacterial bligtit resistance in 
Samba Mahsuri, an elite indica rice variety. Euphytica 160: 411-422. 
doi: 1 0.1 007/s1 0681 -007-9564-6. 

7. Reddy A, Srinivasan T, John V (1980). Breeding for bacterial leaf blight 
resistance in rice. Coimbatore, India: Tamil Nadu Agricultural 
University, pp. 5-13. 

8. Seshu DV (1989). Salient findings from the multllocation evaluation of 
the international rice bacterial blight nursery. Manila, Philippines: 
International Rice Research Institute, pp. 167 - 176. 

9. Leach JE, Rhoads ML, Cruz CV, White FF, Mew TW et al. (1992) 
Assessment of genetic diversity and population structure of 
Xanthomonas oryzae pv. oryzae with a repetitive DNA element. AppI 
Environ Microbiol 58: 2188-2195. PubMed: 1353345. 

10. Nelson RJ, Baraoidan MR, Cruz CMV, Yap IV, Leach JE et al. (1994) 
Relationship between phylogeny and pathotype for the bacterial blight 
pathogen of rice. AppI Environ Microbiol 60: 3275-3283. PubMed: 
16349380. 

11. Leach JE, Leung H, Nelson RJ, Mew T-wW (1995) Population biology 
of Xanthomonas oryzae pv. oryzae and approaches to Its control. 
Current Opinion in Biotechnology 6: 298-304. doi: 
10.101 6/0958-1 669(95)80051 ^. 

12. Adhikari TB, Cruz C, Zhang Q, Nelson RJ, Skinner DZ et al. (1995) 
Genetic Diversity of Xanthomonas oryzae pv. oryzae in Asia. AppI 
Environ Microbiol 61: 966-971. PubMed: 16534980. 

13. Yashitola J, Krishnaveni D, Reddy AP, Sonti RV (1997) Genetic 
diversity within the population of Xanthomonas oryzae pv. oryzae in 
India. Phytopathology 87: 760-765. doi:10.1094/PHYTO.1997.87.7.760. 
PubMed: 18945099. 

14. Shanti ML, George M, Cruz CV, Bernardo M, Nelson R et al. (2001) 
Identification of resistance genes effective against rice bacterial blight 
pathogen in eastern India. Plant Disease 85: 506-512. doi:10.1094/ 
PDIS.2001. 85.5.506. 

15. Lore JS, Vikal Y, Hunjan MS, Goel RK, Bharaj TS et al. (2011) 
Genotypic and pathotypic diversity of Xanthomonas oryzae pv. oryzae, 
the cause of bacterial blight of rice in Punjab State of India. Journal of 
Phytopathology 159: 479-487. doi:10.1 1 1 1/j. 1439-0434.201 1 .01 789.x. 

16. Singh S, Sidhu J, Huang N, Vikal Y, Li Z et al. (2001) Pyramiding three 
bacterial blight resistance genes (xa5, xa13 and Xa21) using marker- 
assisted selection Into indica rice cultivar PR106. Theoretical and 
Applied Genetics 102: 1011-1015. doi:10.1007/s001220000495. 

17. Kauffman H, Reddy A, Hsieh S, Merca S (1973) An improved technique 
for evaluating resistance of rice varieties to Xanthomonas oryzae. Plant 
Dis Rep 57: 537-541. 

18. Sambrook J, Fritsch EF, Maniatis T (1989) Molecular cloning: a 
laboratory manual. Cold Spring Harbor Laboratory. 

19. HampI V, Pavlfcek A, Flegr J (2001) Construction and bootstrap 
analysis of DNA fingerprinting-based phylogenetic trees with the 
freeware program FreeTree: application to trichomonad parasites. Int J 
Syst Evol Microbiol 51: 731-735. doi:10.1099/00207713-51-3-731. 
PubMed: 11411692. 

20. Page RD (1996) Tree View: An application to display phylogenetic 
trees on personal computers. Comput AppI Biosci 12: 357-358. 
PubMed: 8902363. 



21. Salipante SJ, Hall BG (201 1) Inadequacies of minimum spanning trees 
in molecular epidemiology. J Clin Microbiol 49: 3568-3575. dol:10.1128/ 
JCM.00919-11. PubMed: 21849692. 

22. Ellson J, Gansner ER, Koutsofios E, North SC, Woodhull G (2004) 
Graphviz and dynagraph — static and dynamic graph drawing tools. 
Graph drawing software. Springer, pp. 127-148. 

23. Li Z-K, Sanchez A, Angeles E, Singh S, Domingo J et al. (2001 ) Are the 
dominant and recessive plant disease resistance genes similar?: A 
case study of rice R genes and Xanthomonas oryzae pv. oryzae races. 
Genetics 159: 757-765. PubMed: 11606550. 

24. Iyer AS, McCouch SR (2004) The rice bacterial blight resistance gene 
xa5 encodes a novel form of disease resistance. Mol Plant Microbe 
Interact 17: 1348-1354. dol:10.1094/MPMI.2004.17.12.1348. PubMed: 
15597740. 

25. Sugio A, Yang B, Zhu T, White FF (2007) Two type III effector genes of 
Xanthomonas oryzae pv. oryzae control the induction of the host genes 
GsTFIIAyl and OsTFXI during bacterial blight of rice. Proc Natl Acad 
Sci U S A 104: 10720-10725. doi:10.1073/pnas.0701742104. 

26. Boch J, Bonas U (2010) Xanthomonas AvrBs3 family-type III effectors: 
discovery and function. Annu Rev Phytopathol 48: 419-436. doi: 
10.1146/annurev-phyto-080508-081936. PubMed: 19400638. 

27. Chu Z, Yuan M, Yao J, Ge X, Yuan B et al. (2006) Promoter mutations 
of an essential gene for pollen development result in disease resistance 
in rice. Genes Dev 20: 1250-1255. dol:10.1101/gad.1416306. PubMed: 
16648463. 

28. Yang B, Sugio A, White FF (2006) Os8N3 is a host disease- 
susceptibility gene for bacterial blight of rice. Proc Natl Acad Sci U S A 
103: 10503-10508. doi:10.1073/pnas.0604088103. PubMed: 
16798873. 

29. Chen L-Q, Qu X-Q, Hou B-H, Sosso D, Osorio S et al. (2012) Sucrose 
efflux mediated by SWEET proteins as a key step for phloem transport. 

Science 335: 207-211. doi:10.1126/science.1213351. PubMed: 
22157085. 

30. Antony G, Zhou J, Huang S, Li T, Liu B, et al. (2010) Rice xa13 
recessive resistance to bacterial blight is defeated by induction of the 
disease susceptibility gene Os-11N3. The Plant Cell Online 22: 
3864-3876. 

31. Rajeshwari R, Sonti RV (2000) Stationary-phase variation due to 
transposition of novel insertion elements in Xanthomonas oryzae pv. 
oryzae. J Bacterlol 182: 4797-4802. dol:10.1128/JB. 
182.17.4797-4802.2000. PubMed: 10940020. 

32. Yashitola J, Reddy A, Sonti RV (2000) A widely distributed lineage of 
Xanthomonas oryzae pv. oryzae in India may have come from native 
wild rice. Plant Disease 84: 465-469. doi:10.1094/PDIS.2000.84.4.465. 

33. Jeung JU, Heu SG, Shin MS, Vera Cruz CM, Jena KK (2006) 
Dynamics of Xanthomonas oryzae pv. oryzae populations in Korea and 
their relationship to known bacterial blight resistance genes. 
Phytopathology 96: 867-875. doi:10.1094/PHYTO-96-0867. PubMed: 
18943752. 

34. Sundaram RM, Vishnupriya MR, Laha GS, Rani NS, Rao PS et al. 
(2009) Introduction of bacterial blight resistance into Triguna, a high 
yielding, mid-early duration rice variety. Blotechnol J 4: 400-407. doi: 
10.1002/blot.200800310. PubMed: 19253322. 



PLOS ONE I www.plosone.org 



11 



November 2013 | Volume 8 | Issue 1 1 | e81996 



